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ABSTRAGT 

We present a Bayesian framework to account for the magnification bias from both strong and 
weak gravitational lensing in estimates of high-redshift galaxy luminosity functions. We illustrate our 
method by estimating the z ~ 8 UV luminosity function using a sample of 97 Y-band dropouts (Lyman 
break galaxies) found in the Brightest of Reionizing Galaxies (BoRG) survey and from the literature. 

We find the luminosity function is well described by a Schechter function with characteristic magnitude 
of M* = —19.85to;35, faint-end slope of a = —1.72to;29> and number density of log^p 'L*[Mpc“^] = 

—3.00^0 3 ^. These parameters are consistent within the uncertainties with those inferred from the 
same sample without accounting for the magnification bias, demonstrating that the effect is small 
for current surveys at z ^ 8, and cannot account for the apparent overdensity of bright galaxies 
compared to a Schechter function found recently by Bowler et al. (2014a,b) and Finkelstein et al. 

(2014). We estimate that the probability of finding a strongly lensed z ~ 8 source in our sample is 
in the range ~ 3 — 15% depending on limiting magnitude. We identify one strongly-lensed candidate 
and three cases of intermediate lensing in BoRG (estimated magnification fj, > 1.4) in addition to the 
previously known candidate group-scale strong lens. Using a range of theoretical luminosity functions 
we conclude that magnification bias will dominate wide field surveys - such as those planned for the 
Euclid and WFIRST missions - especially at z > 10. Magnification bias will need to be accounted 
for in order to derive accurate estimates of high-redshift luminosity functions in these surveys and to 
distinguish between galaxy formation models. 

Subject headings: galaxies: high-redshift — galaxies: LF — gravitational lensing: strong 


1. INTRODUCTION 

Accurate measurements of the rest-frame UV luminos¬ 
ity function (LF) are crucial for studying the evolution of 
galaxies at high redshift and reconstructing the physics 
and timeline of cosmic reionization. In recent years, sig¬ 
nificant progress has been achieved in measuring the LF 
out to z ~ 8 and beyond based on images taken with 
the Hubble Space Telescope in the deep legacy fields, the 
Hubble Frontier Fields and through parallel programs 
(e.g., Bouwens et al. 2007, 2011; Bradley et al. 2012; 
Oesch et al. 2012; Robertson et al. 2013; McLure et al. 
2013; Schenker et al. 2013; Schmidt et al. 2014a; Bradley 
et al. 2014; Bouwens et al. 2014; Finkelstein et al. 2014; 
Zitrin et al. 2015; Goe et al. 2015). 

From many of these surveys it appears the LF at z < 6 
is well fit by a Schechter (1976) function with a power- 
law slope at faint luminosities and an exponential drop at 
the bright end, where it is expected that feedback reduces 
star-formation in the most massive galaxies (Somerville 
et al. 2012) and dust extinction may reduce the UV flux 
of galaxies (Cai et al. 2014). The evolution of the LF is 
expected to be driven by these processes and the evolu¬ 
tion of the underlying halo mass function. It is so far un- 
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established which processes dominate the evolution and 
whether there are signification changes in the physical 
conditions of galaxies forming at high redshifts. 

Recent studies by Bowler et al. (2014a,b) and Finkel¬ 
stein et al. (2014) claimed an over-abundance of galaxies 
at the bright end of the z > 6 LF when compared to 
the fit of a Schechter (1976) function, although Bouwens 
et al. (2014) found no evidence for a departure from 
a Schechter-like form at z ~ 4 — 8, largely analyzing 
the same data sets. An over-abundance of bright galax¬ 
ies may also be apparent in smaller surveys (e.g., Ono 
et al. 2012; Hathi et al. 2012; Finkelstein et al. 2013). 
If the departure from an exponential cutoff is confirmed 
by future observations, this may be an indication of the 
changing astrophysical conditions of high-redshift galax¬ 
ies. However, another possible explanation is that the LF 
remains intrinsically with a Schechter form and the over¬ 
abundance of bright galaxies is caused by gravitational 
lensing magnification bias, which has been predicted to 
be significant for galaxies at z > 8 (Wyithe et al. 2011). 

While it has long been recognized that the gravita¬ 
tional lensing effect can be exploited in order to probe 
intrinsically faint galaxies - in particular behind massive 
clusters of galaxies at moderate redshift - (e.g., Franx 
et al. 1997; Ellis et al. 2001; Schmidt et al. 2014b; Bowler 
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et al. 2014b; Zitrin et al. 2015; Atek et al. 2015; Coe et al. 
2015), the effect in blank fields is much less well appre¬ 
ciated. 

In fact, gravitational leasing affects all lines of sight, 
as the trajectory of every photon in the universe is per¬ 
turbed by the inhomogeneous foreground mass distribu¬ 
tion. Though the effect is generally not as strong as in 
the fields of massive clusters of galaxies, even so-called 
blank field surveys are affected by gravitational leasing 
(weak, intermediate, or strong). In practice, owing to the 
leasing effect, flux-limited surveys include sources that 
should be below the sample threshold, but have been 
magnified into the sample. Furthermore, gravitational 
leasing changes the relation between observed solid an¬ 
gle and cosmic volume with respect to that expected for a 
perfectly homogeneous universe. At fixed detector field- 
of-view the intrinsic solid angle observed is smaller for 
magnification ^ > 1 and vice versa. This phenomenon is 
called magnification bias (e.g.. Turner et al. 1984; Wyithe 
et al. 2001, 2011) and it can change the shape of the ob¬ 
served LF. Thus, it needs to be accounted for in order to 
derive accurate intrinsic LFs from flux-limited samples. 

The main aim of this paper is to improve the estima¬ 
tion of the intrinsic UV LF at high redshift by developing 
a formalism to take into account the magnification bias. 
Our new formalism improves on previous work in sev¬ 
eral ways: we extend the analytic strong lensing model 
of Wyithe et al. (2011) to include the redshift evolution 
of the deflector population, and we develop a technique 
to treat the intermediate lensing regime and introduce 
a framework to include weak lensing effects, neither of 
which have been systematically accounted for in any pre¬ 
vious estimates of the LF. Furthermore, by providing 
probability distribution functions for the magnification 
of each dropout and empty field, our formalism can be 
directly included in any Bayesian LF parameter estima¬ 
tion, thus allowing for a rigorous derivation of the related 
uncertainties. 

We present two applications of our formalism. The first 
application is the interpretation of the z ~ 8 dropouts 
found by the Brightest of Reionizing Galaxies Survey^, 
(hereafter BoRG, Trenti et al. 2011). After estimating 
the fraction of sources in BoRG that are multiply imaged 
and presenting one strongly-lensed candidate and three 
candidate systems with magnification fi > 1.4, we use the 
extended sample presented by Schmidt et al. (2014a) to 
derive the LF including the effects of magnification bias. 
For this we extend the Bayesian formalism introduced 
by Schmidt et al. (2014a) by including a term describing 
the likelihood for magnification of high-redshift sources 
for each field, and marginalize over the range of possible 
magnifications. 

The second application of our formalism is a set of pre¬ 
dictions for the modification of the LF at 8 < z < 16, 
where JWST will detect dropouts (Gardner et al. 2006), 
by using a variety of possible LFs based on theoreti¬ 
cal models (Munoz 2012; Behroozi & Silk 2015) and ex¬ 
trapolations of lower redshift data (Bouwens et al. 2014; 
Finkelstein et al. 2014). With our formalism we can give 
a quantitative assessment of how magnification bias will 
affect future surveys. 

The paper is organized as follows. In Section 2 we 
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briefly describe the BoRG survey and the data used in 
this paper. In Section 3 we introduce the relevant the¬ 
oretical background for gravitational lensing and mag¬ 
nification bias. In Section 4 we develop a semi-analytic 
framework, based on that in Wyithe et al. (2011) to study 
the magnification bias due to strong and intermediate 
gravitational lensing. In Section 5 we use the reconstruc¬ 
tion of lines-of-sight in cosmological simulation data to 
investigate weak lensing. The Bayesian inference for the 
determination of the intrinsic LF is introduced in Sec¬ 
tion 6 and presented in more detail in Appendix A. The 
results are presented and discussed in Section 7. A brief 
summary is given in Section 8. 

All magnitudes are AB magnitudes and a standard 
concordance cosmology with flm = 0.3, IIa = 0.7, and 
h = 0.7 is assumed. The Millennium Simulation uses a 
cosmology with = 0.25, = 0.75, and h = 0.73, 

which is used to estimate the weak lensing magnification. 
We assume the difference between these two cosmologies 
is negligible for our purposes. 

2. DATA 

This paper estimates the z ~ 8 LF using 38 bright 
Lyman Break galaxies selected from the BoRG survey 
and 59 fainter dropouts taken from deep legacy fields 
(in HUDF09 and the WFC3/IR wide area Early Re¬ 
lease Science). The BoRG survey is described briefly 
in Section 2.1, but we refer to Trenti et al. (2011, 2012); 
Bradley et al. (2012) and Schmidt et al. (2014a) for fur¬ 
ther details. The deep legacy data are described by 
Bouwens et al. (2011). Additionally, we used data of 
galaxies with spectroscopically determined velocity dis¬ 
persions to estimate the velocity dispersion of the fore¬ 
ground BoRG galaxies (described in Section 2.2). In Sec¬ 
tion 2.3 we give an overview of the simulated data used 
in the analysis of weak lensing. 

2.1. The BoRG Survey 

The ongoing BoRG survey is a pure-parallel imaging 
program with the HST WFC3. The current survey covers 
^ 350 arcmin^ divided into 71 independent fields located 
randomly on the sky. This reduces cosmic variance be¬ 
low the level of statistical noise (Trenti & Stiavelli 2008; 
Bradley et al. 2012). The photometry is in the visual 
and near-infrared, primarily using the four HST WFC3 
filters F606W, F098M, F125W, and F160W (commonly 
referred to as V-, Y-, J-, and H-bands respectively). The 
z ^ 8 BoRG survey consisted mainly of HST programs 
GO/PAR 11700 and GO/PAR 12572 (PI: Trenti) and 
includes a small additional number of coordinated par¬ 
allels from GOS-GTO. 53 core BoRG fields are comple¬ 
mented by other archival data including 8 fields from 
GO/PAR 11702 (PI: Yan, Yan et al. 2011) and 10 GOS- 
GTO fields, which used the F600LP-band instead of the 
F606W-band. The BoRG survey is the largest current 
survey of Y-band dropouts by solid angle. 

The z ^ 8 galaxy candidates were identified from 
Y-band dropouts, full details of the selection criteria 
used to find dropouts are described in Schmidt et al. 
(2014a). The BoRG survey detected 38 Lyman break 
galaxy (LBG) candidates at z ~ 8 with S/N > 5 in 
the J-band, of which 10 have S/N > 8 (Bradley et al. 
2012; Schmidt et al. 2014a). We use the 5cr sample of 
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objects in this work. Throughout this work we will as¬ 
sume 42% of the selected BoRG dropouts are contam¬ 
inants (usually z ~ 2 interlopers, see e.g., Hayes et al. 
2012; Bouwens et al. 2013). This is the fiducial contam¬ 
ination fraction for the BoRG sample and was shown to 
be robust in the estimation of the LF by Bradley et al. 
(2012); Schmidt et al. (2014a). By definition we can¬ 
not determine which specific sources are contaminants 
without further photometry and spectroscopy, but our 
rigorous Bayesian method to determine the LF allows 
us to accurately estimate the LF parameters accounting 
for the presence of random contaminants (Schmidt et al. 
2014a). 


0 ~ 12. While the Millennium Simulation contains ha¬ 
los from very high redshift, it has a box length of only 
500 Mpc h~^. The comoving distance in the universe 
to z = 1 is 2390 Mpc so it is necessary to build 
lightcones with the galaxies correctly distributed in co¬ 
moving volumes (see Blaizot et al. (2005) and Kitzbichler 
& White (2007) for a thorough discussion of generating 
mock lightcones). 

These lightcones were generated using the semi- 
analytical galaxy formation model of Guo et al. (2011), 
and photometric properties were calculated using the 
stellar population synthesis code by Maraston (2005) 
which can be applied at high redshift. 


2.2. Massive Foreground Galaxies Acting as Deflectors 

In Section 4.2 we estimate the velocity dispersions of 
strong lens candidates in the BoRG fields by comparing 
their photometry with similar early-type galaxies which 
have both HST photometry and spectroscopically deter¬ 
mined velocity dispersions (Treu et al. 2005; Belli et al. 
2014a,b). We divided the galaxy samples into three large 
redshift bins in order to account for the position of the 
4000A break in the filters at higher redshifts. 

In the range z < 1 we used a sample of 165 spheroidal 
galaxies from Treu et al. (2005) with photometry from 
the Great Observatories Origins Deep Survey North 
(GOODS-N, Bundy et al. 2005). For z > 1 we use a sam¬ 
ple of 66 massive quiescent galaxies, presented by Belli 
et al. (2014a,b), which were selected from HST photo¬ 
metric catalogs of objects in the GOSMOS, GOODS and 
Extended Groth Strip (EGS) fields (Grogin et al. 2011; 
Koekemoer et al. 2011; Windhorst et al. 2011). We used 
an aperture correction to rescale observed velocity dis¬ 
persions, (Tobs) to (Je, the velocity dispersion within one 
effective radius. Re- 

We follow Belli et al. (2014a) and used the model of 
van de Sande et al. (2013) which proposes a constant 
rescaling: 

(Te = l.OhCTobs (I) 

For galaxies at z < 1 (the Treu et al. (2005) sample), we 
used the model of Cappellari et al. (2006): 

_ 

— \ J ^obs \^) 

where the slit size, R is the 1" aperture on Keck DEIMOS 
(Treu et al. 2005). 

The reference photometry used for the individual sam¬ 
ples differ. As listed in Table 1 we use HST F606W for 
galaxies at z < 0.5, HST F850LP from Treu et al. (2005) 
(converted to F098M through linear interpolation) for 
galaxies at 0.5 < z < 1.0, and HST F160W for galaxies 
at z > 1. 


2.3. The Millennium Simulation 

In Section 5 we describe our method to generate weak 
leasing probability density functions (PDFs) by recon¬ 
structing simulation data along the line-of-sight to z ~ 8. 
Due to the very high redshift of our sources, it was neces¬ 
sary to use simulation data containing halos out to red- 
shifts above 5. 

We used 24 1.4 x 1.4 square degree simulated lightcones 
built by Henriques et al. (2012) from the Millennium Sim¬ 
ulation (Springel et al. 2005) which contain halos out to 


3. THEORETICAL BACKGROUND 

In this section we summarize the relevant theory for 
the galaxy LF, strong and weak gravitational leasing, 
and magnification bias. 

3.1. Galaxy Luminosity Function 

When a simply parametrized form is needed, we de¬ 
scribe the LF by a Schechter function (Schechter 1976): 



where L* marks the characteristic break in the LF, 4'* 
is the characteristic density at that luminosity and a is 
the power-law exponent slope of the faint end. 


3.2. Strong Sensing 

If the line-of-sight to a background source is closely 
aligned with a massive foreground object, e.g. a cluster 
or single massive galaxy, gravitational lensing can pro¬ 
duce multiple observed images of the source (Schneider 
et al. 1992, 2006). Multiple imaging signifies the regime 
of strong gravitational lensing lensing. 


3.2.1. Singular Isothermal Sphere 

Strong gravitational lenses are commonly modeled as 
Singular Isothermal Spheres (SIS), which provides a con¬ 
venient analytic form to describe the mass profiles of 
massive galaxies (e.g. Treu 2010, and references therein). 
The scale of image separation is characterized by the Ein¬ 
stein radius of the lens: 

^'er(ct,z) = 47r^ (4) 

where Dig and Dg are the angular diameter distances be¬ 
tween the lens and source, and from the observer to the 
source respectively, a is the velocity dispersion of the 
lens galaxy, and c is the speed of light. Velocity disper¬ 
sion is the most important property for determining the 
strength of a strong gravitational lens as it scales with 
the mass of the dark matter in the system (Turner et al. 
1984; Schneider et al. 2006; Treu 2010). 

The magnification, yt, due to an SIS lens is given by: 


S = 


\ 0 \ 

\d\ — ^ER 


( 5 ) 


where 0 is the distance between the lens and the source 
in the image plane. An SIS lens can produce two images, 
with the brighter one having magnification /i > 2, or one 
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image with magnification fx < 2. The case of multiple 
imaging is referred to here as strong lensing. In this pa¬ 
per we refer to images with 1.4 < /i < 2 as intermediate 
lensing. 

3.2.2. Multiple Image Optical Depth 

The optical depth is the cross-section for a galaxy at 
redshift zs to be multiply imaged (i.e. strongly lensed) 
by a foreground galaxy ai z^: it is the fraction of the 
sky covered by the Einstein radii of all intervening de¬ 
flectors at redshifts zl- Following standard practice and 
assuming SIS deflectors, Wyithe et al. (2011) defines it 
as: 


r^s r 

= / dzL / da ^{a,ZL) {I+ZlT 

( 6 ) 

where ^{a, zl) is the velocity dispersion function of the 
deflectors, Dp is the angular diameter distance to zl, 
and t is time. Without the magnification bias, the optical 
depth gives the probability of a high-redshift source being 
multiply imaged. 

3.3. Weak Lensing 

Weak gravitational lensing is the deflection of light that 
causes the magnification and distortion of an observed 
source, but without producing multiple images. There 
are no empty lines-of-sight in the universe, so all light 
traveling to us has been deflected some amount by inter¬ 
vening mass (Hilbert et al. 2007). Whilst it is impossi¬ 
ble to determine the exact effect on individual observed 
sources, it can be done in a statistical sense and is impor¬ 
tant to quantify this effect for our high-redshift sources. 

The lens equation can be constructed for an arbitrary 
number of lens planes due to an ensemble of deflectors 
along the line-of-sight (Hilbert et al. 2009; McCully et al. 
2014). The magnification of a source in a multiplane 
system is a function of the total convergence and total 
shear experienced. Hilbert et al. (2009) showed to first 
order that the total convergence and shear are the sum of 
the individual contributions from each object along the 
line-of-sight: 

p ( 7 ) 

(1 - \ 

The convergence, Ki, and shear, 7^, of each object are 
determined by the lens model. 

3.4. Magnification Bias 

The gravitational lensing of a source with luminosity 
T in a solid angle Vl of sky has two effects. The observed 
luminosity is magnified by a factor /x and sources are now 
distributed over a magnified solid angle /ill. In a flux- 
limited sample intrinsically low luminosity sources can 
be magnified above the survey limit, while the number 
density of sources can decrease for a given observed solid 
angle. 

Since the faint end of the LF of high-redshift LEG 
galaxies is so steep, in regions around large low-redshift 
deflectors we may observe an excess of intrinsically faint 
high-redshift sources. These effects are known as the 
magnification bias and will effect our inferences about 
the population and LF of high-redshift galaxies. 


If it were possible to observe all galaxies in the universe 
without the magnification bias the probability of a high- 
redshift galaxy being strongly lensed is purely given by 
the optical depth, Tm (Section 3.2.2). However, magni¬ 
fication of more numerous intrinsically faint sources into 
our surveys implies that we do not observe the true pop¬ 
ulation of galaxies with luminosity. The magnification 
bias increases the probability that a sample of observed 
high-redshift sources have been gravitationally lensed. 

The magnification bias for sources with observed lu¬ 
minosities above Tiim in a flux-limited sample is given 
by: 


/;::‘ibp(b)«(>^) 

N{> Lliny) 


( 8 ) 


assuming that each source could be magnified between 
Mmin and /imax- Where p(/i) is the probability distribu¬ 
tion for magnification of a source and N{> Liim) is the 
integrated galaxy LF (Wyithe et al. 2011). 

The true probability of a high-redshift source being 
multiply imaged is BTm- Therefore, using B it is possible 
to find the fraction of galaxies at a given redshift in a 
flux-limited sample that are multiply-imaged: 

i^ynult = „ -r (9) 

BTjn + B'{1 - Tm) 


We assume that B', the bias for galaxies to not be mul¬ 
tiply imaged is close to unity. 

If the survey limit is brighter than the characteristic 
apparent magnitude of the observed sample the magnifi¬ 
cation bias is expected to be large, as a large fraction of 
the observed sources are likely to be intrinsically fainter 
sources magnified above the detection threshold of the 
survey. 

We can compute the gravitationally lensed LF, includ¬ 
ing strong and weak gravitational lensing: 


«'mod(F) = (l-Ty„)— ^ ) 

Mdemag \ Mdemag / 


+ Tr 


, r dfx(-) ( 10 ) 

Jo L \L J 


Where /Xdemag < 1 is introduced such that the mean mag¬ 
nification over the entire sky is unity (Pei 1995; Wyithe 
et al. 2011) and p(/x) is the full probability density for 
magnification of a high-redshift source, as above. For a 
Schechter LF, the gravitationally lensed LF is predicted 
to exhibit a ‘kick’ in the bright end (e.g., Wyithe et al. 
2011) due to a pile-up of brightened galaxies, whereas 
at the faint end the magnification of flux is balanced by 
the loss of number density (for faint-end slope a ~ —2, 
Blandford & Narayan 1992) so there is no distortion, even 
if many strongly lensed faint sources are observed. 


4. STRONG AND INTERMEDIATE LENSING 

In this section we compute the probability that the 
2; ~ 8 dropouts are affected by strong and intermedi¬ 
ate lensing. First, in Section 4.1 we compute the strong 
lensing optical depth and the probability that a z ~ 8 
source is multiply imaged by foreground massive ellip¬ 
tical galaxy deflectors. We account for evolution of 
the deflector population based on the observed stellar 
mass function. In Section 4.2 we describe our method 



LF Magnification Bias 


5 


to identify sources in the intermediate leasing regime 
(1.4 < <2). In order to identify these sources, we esti¬ 

mate the lensing strength of massive foreground galaxies 
based on HST photometry and an empirical calibration 
of the Faber & Jackson (1976) relation. A candidate 
strongly leased dropout in the BoRG fields was presented 
by Barone-Nugent et al. (2013), in this paper one more 
candidate multiply-imaged dropout (/x > 2) is found, and 
three dropouts may experience significant intermediate 
magnification. We detail their properties in Table 2. 

4.1. Strong Lensing by an Evolving Deflector Population 

In order to compute the strong lensing optical depth 
and multiple image probability, we follow Wyithe et al. 
(2011) and use a simple SIS lensing model (see Sec¬ 
tion 3.2.1) with a flat cosmology. Strong lenses are as¬ 
sumed to be uniformly distributed in the universe and 
we can calculate the probability of encountering a strong 
lens along the line-of-sight to a high-redshift source, i.e. 
the lensing optical depth (see Section 3.2.2). By consid¬ 
ering the number of galaxies observed above a certain 
flux limit we can calculate the magnification bias factor, 
B, from Equation (8), assuming a Schechter luminos¬ 
ity function (Equation (3)). For these calculations we 
use the z ~ 8 LF inferred by Schmidt et al. (2014a), 
with a characteristic magnitude of M* = —20.15^g'3g, 
faint-end slope of a = —1.87 ^q 26, number density 

of logiQ 4'*[Mpc“^] = —3.24^0^4- We marginalize over 
the entire MCMC chain for each of the Schechter param¬ 
eters. 

In their calculation of the optical depth Wyithe et al. 
(2011) used the local velocity dispersion function as mea¬ 
sured by SDSS (Choi et al. 2007). As most strong lenses 
occur at z < 1.5 (Fassnacht et al. 2004; Treu 2010), 
Wyithe et al. (2011) assumed that the velocity dispersion 
function does not evolve with redshift for massive galax¬ 
ies. This is consistent with studies of the velocity dis¬ 
persion function out to z ~ 1 (e.g., Chae 2010; Bezanson 
et al. 2012). However, significant galaxy growth and evo¬ 
lution is observed from z > 1 as structure forms (van de 
Sande et al. 2013; Belli et al. 2014a), and we can improve 
the accuracy of the model by allowing the parameters of 
the velocity dispersion function for massive ellipticals to 
evolve with redshift. Introducing redshift evolution is 
expected to reduce the optical depth (Barkana & Loeb 
2000 ). 

The dashed blue line in the left panel of Figure 3 shows 
the probability that the source has been multiply imaged 
as a function of lens redshift for a source at z ~ 8, cal¬ 
culated using Equation (6). The distribution is strongly 
peaked at zz, ~ 1, but there is a significant probability 
that zl > 1.5. Only 48% of the contribution to the op¬ 
tical depth for strong lensing occurs at zl < 1-5. We 
find that 90% of lensing occurs within a lens redshift of 
Zl < 3.5. Therefore, in order to account for most of the 
optical depth we need to find the form of the velocity dis¬ 
persion function out to z = 3 ~ 4 where the galaxy popu¬ 
lation is significantly different from recent times (Bundy 
et al. 2005; Muzzin et al. 2013; van de Sande et al. 2014). 

Several studies have investigated the evolution of the 
velocity dispersion function out to z ~ 1.5 (e.g., Chae 
2010; Bezanson et al. 2011, 2012, 2013). These works 
are consistent with no evolution, but have large uncer¬ 


tainties. Measurements of velocity dispersion beyond 
z > 2 are very difficult as the brightest emission lines fall 
within near-IR atmospheric absorption regions (Kriek 
et al. 2006; Belli et al. 2014b). 

Therefore, we estimate the evolution of the velocity dis¬ 
persion function at high redshift based on the evolution 
of the stellar mass function, a related quantity that has 
been well-measured at z > 2. We convert the stellar mass 
function into the velocity dispersion function by means of 
the well-known correlation between stellar velocity dis¬ 
persion (cr) and stellar mass (Mgteu) taken from Auger 
et al. (2010): log(cr[km s“^]) = pM — lip -|- q, where 
p = 0.24±0.02, q = 2.34±0.01 and M = log (M^teii/Mo). 
This relation was derived for massive lens galaxies with 
high velocity dispersions, which will be the strongest con¬ 
tribution to the optical depth as r ~ ct'^. 

High-redshift galaxies are observed to have higher ve¬ 
locity dispersions at fixed mass than in the local uni¬ 
verse (e.g., van de Sande et al. 2013; Belli et al. 2014a; 
Bezanson et al. 2015). Thus the stellar mass-velocity 
dispersion relation is expected to evolve with redshift. 
Following van de Sande et al. (2013) we expect evolu¬ 
tion of the form (ct/cto) oc (1 -I- z)^, where (Jq is the ex¬ 
pected velocity dispersion at z ~ 0. In Figure 1 we plot 
publicly available data from van der Wei et al. (2008); 
van Dokkum et al. (2009); Newman et al. (2010); Toft 
et al. (2012); Bezanson et al. (2013); van de Sande et al. 
(2013); Belli et al. (2014a,b) and fit a relation of this 
form for all galaxies with estimated stellar masses be¬ 
tween 10.8 < log(Msteii/A70) < 12.0, and measured ve¬ 
locity dispersion a > 200 km s“^ as this was the region 
where the Auger et al. (2010) relation was derived. We 
find fl = 0.20 ± 0.07. Our result is lower than the re¬ 
sult from van de Sande et al. (2013) because we use the 
Auger et al. (2010) stellar mass-velocity dispersion rela¬ 
tion for massive lens galaxies as CTq, whereas van de Sande 
et al. (2013) compare to a dynamical mass-velocity dis¬ 
persion relation. As demonstrated in van de Sande et al. 
(2013) Msteii/Afdyn increases with redshift, so will reduce 
the evolution we find compared to that in van de Sande 
et al. (2013). If we consider the same galaxy sample and 
fit both our relation derived from stellar masses and the 
van de Sande et al. (2013) dynamical mass relation, and 
include the evolution in Msten/A7dyn7 our results are con¬ 
sistent. We note that because the optical depth depends 
on velocity dispersion to the fourth power, the form of 
the velocity dispersion function at z > 2 is the greatest 
source of uncertainty in the calculation of optical depths. 

The stellar mass function can be described by a 
Schechter function (e.g., Muzzin et al. 2013): 


$g(M) = (In 10) io(“-“s)(i+“s) exp 


-10 


M-Ml 


The characteristic stellar mass is given by Mg = 
log (MJ^gjj/M©), <l>g is the characteristic density normal¬ 
ization, and as is the low-mass-end slope. 

In order to model the redshift evolution of the stellar 
mass function, we use publicly available data on quies¬ 
cent galaxies at z < 4 from the COSMOS/UltraVISTA 
Survey (Muzzin et al. 2013). They derive the best- 
fit single Schechter function parameters for the stellar 
mass function as a function of redshift. Their stellar 
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Fig. 1.— Redshift evolution of massive galaxy velocity dispersion, 
relative to the velocity dispersion estimated from inferred stellar 
masses via the Auger et al. (2010) relation. We find evolution of the 
form ((t/(to) oc (1 + where (tq is the velocity dispersion 

estimated using the stellar mass-velocity dispersion relation from 
Auger et al. (2010). We plot the mean linear fit (black line) and 
the la confidence region (gray shaded region). 

mass function parameters for quiescent galaxies, allow¬ 
ing for evolution of as, are plotted as a function of red- 
shift in Figure 2. We assumed the redshift evolution 
X = Xo(l + z)°‘, where X represents the stellar mass 
function Schechter parameters and Xq represents the val¬ 
ues a,t z = 0. 

We used a Bayesian MCMC linear fitting method to 
fit this functional form to the data, and plot the mean 
and one standard deviation confidence fits in Figure 2. 
There is significant evolution in However, there is 
also large uncertainty in the evolution of due to the 
spread of the data. We ignore evolution in the low-mass- 
end slope, since the lensing effect is dominated by the 
most massive galaxies. We also ignore evolution in Mg, 
for which the evolution appears non-negligible but it has 
little effect on Equation (11). The redshift-dependent 
velocity dispersion function obtained in this way becomes 


^{a,z) =p 


cr(l -I- z)^ 


/ a \P Ai+“s) 

f (T \P ^ 



\a / 

\a / 


( 12 ) 

with p = 0.24 ± 0.02, /3 = 0.20 ± 0.07 (obtained 
from the evolution of velocity dispersion in Figure 1), 
$*(z) = 3.75 ± 2.99 X 10-3(1 -b 2)-2-46±o.53 Mpc-^, 
as = -0.54 ± 0.32 and cr* = 216 ± 18 km s"!. This 
was derived using the stellar mass-velocity dispersion re¬ 
lation above (Auger et al. 2010), including the scatter in 
the relation. At z = 0 recent well-measured velocity dis¬ 
persion functions (e.g., Sheth et al. 2003; Choi et al. 2007) 
are within the uncertainties of this redshift-evolving re¬ 
lation, showing that our inferred evolution is consistent 
with direct measurements where they overlap. 

Using this redshift-dependent velocity dispersion func¬ 
tion we compute the optical depth for strong lensing, 
and the distribution of the optical depth with lens red¬ 
shift. In the left panel of Figure 3, now using the red¬ 
shift evolving deflector population from Equation (6) and 
Equation (12), we see that the majority of the contribu- 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 

log(l + z) 


Fig. 2.— Redshift evolution of the best-fit single Schechter func¬ 
tion parameters from Muzzin et ah (2013) for the stellar mass func¬ 
tion of quiescent galaxies, allowing for evolution of a5. Fits of the 
form Xo{l-\- z)°‘ are plotted: the solid lines show the mean fit, dot¬ 
ted lines show the la error on the data. Only shows significant 
evolution with redshift. 


tion to the optical depth is from lens galaxies at z < 1.5, 
which agrees with current observations of lensed high- 
redshift dropouts (Barone-Nugent et al. 2013; Schmidt 
et al. 2014b; Atek et al. 2015). In the right panel of Fig¬ 
ure 3 we plot the optical depth as a function of source 
redshift and find that including the redshift evolution of 
the deflector population reduces the optical depth at high 
redshift compared with the work in Wyithe et al. (2011) 
as expected by theoretical predictions (Barkana & Loeb 
2000), and it appears to start to flatten by zs ~ 10. 
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Fig. 3.— (Left) Contribution to the optical depth for a source at 2 : ~ 8 to be multiply imaged as a function of the lens redshift, (solid 
black line) calculated using Equation (6), including the evolution of the deflector population with redshift (Section 4.1), for comparison we 
plot the contribution for a constant comoving density of lens galaxies (dashed blue line, Wyithe et al. 2011). (Right) Optical depth for 
multiple imaging as a function of source redshift, including evolution of the deflector population (solid black line). The gray shaded regions 
show the Icr uncertainty bounds on the optical depth and its distribution, given the uncertainties in velocity dispersion and stellar mass 
evolution described in the text. The optical depth without redshift evolution of lens galaxies is also plotted for comparison (dashed blue 
line, Wyithe et al. 2011). 


Our estimated optical depth at z < 8 is in good agree¬ 
ment with values derived by an independent method by 
Barone-Nugent et al. (2015), and consistent with Wyithe 
et al. (2011) for z < 8. We note the optical depths 
presented in Barone-Nugent et al. (2015) are marginally 
higher than the results of this paper, but we can recover 
their optical depth using a steeper evolution of cr(z). It 
is clear that the uncertainty in the evolution of velocity 
dispersion, which is the best indicator of the mass of lens 
galaxies, provides the largest uncertainty in determining 
the optical depth. 

Finally, we compute the probability that high-redshift 
galaxies in a flux-limited sample have been multiply im¬ 
aged. This is shown in Figure 4 as a function of limiting 
magnitude for each of the BoRG fields. As expected, the 
probability that a source in each field is multiply imaged, 
Alnuit (Equation (9)) increases with the survey limiting 
magnitude, owing to the magnification bias. We esti¬ 
mate 3-15 % of observed sources brighter than M* have 
been strongly lensed, this is consistent with the results 
of Barone-Nugent et al. (2015) who use an independent 
method to infer the lensed fraction. 

4.2. Identifying Significantly Magnified Sources 

Whilst all the fields are subject to weak leasing, it is 
necessary to establish which of the individual sources ex¬ 
perience multiple-imaging (/i > 2), or are close enough to 
a deflector to experience an intermediate magnification 
(1.4 < yL <2). We expect strong leasing evens to be rare, 
but possible given the size of the BoRG survey. Among 
the BoRG sources, Barone-Nugent et al. (2013) presented 
a candidate strongly-lensed system in borg_0440-5244 
(for naming conventions see Bradley et al. 2012). The 
candidate appears to be lensed by a foreground group 
with an Einstein radius of 1.49", corresponding to a ve¬ 
locity dispersion of ~ 300 km s“^, producing a magni¬ 
fication of 3.7 ± 0.2 of the dropout. In this Section we 
describe a method to identify other potentially lensed 
sources in the catalogs and illustrate how to account for 



Fig. 4.— Multiply-imaged fraction (see Equation (9)) for 2 ~ 8 
sources brighter than the J-band limiting magnitude in each of the 
BoRG fields, as a function of the UV characteristic magnitude, M*, 
including the evolution of the deflector population (Section 4.1). 
The probability of a high-redshift source being multiply imaged 
increases as the survey magnitude limit becomes brighter than M*. 
We expect very few intrinsically bright sources, so any bright source 
has a high likelihood of being significantly magnified according to 
the magnification bias. We have used the full MCMC chain for M* 
from Schmidt et al. (2014a) and plot the mean value with errorbars 
of one standard deviation. The optical depth, the probability of 
multiple imaging without including the magnification bias factor, 
B (Section 3.4), is plotted as the green dashed line. 

them systematically when estimating the LF. 

For computational speed, we considered as potential 
deflectors only z < 3 objects within 18 arcseconds of 
the z ~ 8 dropouts in each field (the typical Einstein 
radius is of order 1-2 arcseconds for massive galaxies). 
The key quantity that we need to estimate the leasing 
strength is the velocity dispersion (Turner et al. 1984; 
Treu 2010). Thus for every galaxy sufficiently close to a 
dropout, we estimate their velocity dispersions by com¬ 
paring their photometry with that of samples of similar 
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objects with spectroscopically-determined velocity dis¬ 
persions. We selected galaxy samples with HST photom¬ 
etry in bands used in BoRG in order to estimate velocity 
dispersion based on our own photometry. As a compar¬ 
ison sample, we used data from Treu et al. (2005) and 
Belli et al. (2014a,b), as described in Section 2.2. 

As described in Section 4.1, the velocity dispersion- 
stellar mass relation is believed to evolve weakly with 
redshift since z ~ 2 (e.g., van de Sande et al. 2013; Belli 
et al. 2014b; Bezanson et al. 2015), and galaxies will be 
intrinsically brighter at higher redshift due to younger 
stellar populations (e.g. Treu et al. 2005). We account 
for this by fitting an evolving Faber & Jackson (1976) 
relation to the comparison sample of the form L oc + 
z)^. 

In practice, we bin the data in redshift, and fit a func¬ 
tion of the form log a = —0.1m -I- a log(l + z) + b using a 
Bayesian MCMC estimation where a is the velocity dis¬ 
persion in km s“^, m is apparent magnitude in a given 
band, z is galaxy redshift, and a and b are constants. 
We restrict our fit to galaxies with a measured veloc¬ 
ity dispersion of at least 200 km s“^, where samples are 
less affected by incompleteness and selection effects. We 
present the estimated parameters in Table 1 and fits to 
the data are shown in Figure 5. 


TABLE 1 

Correlation between velocity dispersion, redshift and 

APPARENT MAGNITUDE 


Redshift 

Band (m) 

a 

b 

z < 0.5 

F606W 

2.26 ±0.79 

4.08 ±0.12 

0.5 < z < 1.0 

F089M 

0.93 ±0.13 

4.20 ±0.03 

z > 1.0 

F160W 

1.02 ±0.15 

4.12 ±0.05 


Note. - Fits of the form logcr = —0.1m -I- a log(l + z ) + b 


The posterior probability distribution function of Ein¬ 
stein radii for each object are found using Equation (4), 
sampling over the full MCMC chain for the velocity dis¬ 
persion. The redshifts of the objects were determined 
using the Bayesian Photometric Redshifts (BPZ) code, 
using a flat prior and the default parameters and tem¬ 
plates (Benitez et al. 2004; Coe et al. 2006). All photo¬ 
metric redshifts for relevant foreground galaxies are well- 
fit by BPZ and have uncertainties in photometric red¬ 
shift < 15%. The PDF for magnification, p(/i) is found 
by computing the magnification, fj, (Equation (5)), at 
the position of the dropout given the distribution of Ein¬ 
stein radii found for each foreground object using the dis¬ 
tribution for its velocity dispersion, Uinf estimated from 
the fits in Table 1. The greatest source of error in this 
procedure is the magnitude-velocity dispersion-redshift 
relation: uncertainties in magnitude and redshift deter¬ 
mination have small effects on the magnification PDFs 
in comparison to the uncertainty in velocity dispersion. 

When the mean magnification produced by such a fore¬ 
ground object exceeds n = lA we use the magnification 
PDF derived from the above procedure and treat the 
dropout as described in Section 6.2.2 in our calculations 
of the LF. 

Using this method, we find one of the dropouts 
(borg_0436-5259_1233, presented in Bradley et al. 2012) 
has a magnification probability distribution consistent 
with strong lensing. This dropout is shown in the top 


left panel of Figure 6 and its estimated lensing proper¬ 
ties are given in Table 2. The dropout appears to be 
magnified by a large galaxy at z ~ 0.40 with estimated 
velocity dispersion 294 ±47 km s“^ (estimated from pho¬ 
tometry via the empirical relation presented in Table 1). 
We estimate its magnification to be /r = 2.05 ± 0.52, 
the large uncertainty is due to the uncertainty in the 
relationship between apparent magnitude and velocity 
dispersion. If the dropout is indeed strongly lensed the 
counter image would be almost directly behind the center 
of the lens galaxy, and will be demagnified according to 
Equation (5), unfortunately making it impossible to de¬ 
tect. The dropout is very faint (mj = 27.0 ±0.2) and no 
significant elongation is detected in any of the observed 
bands but this dropout would be an excellent object for 
further investigation. 

Three of the dropouts (borg_1301±0000_160, 
borg_1408±5503 and borg_2155-4411_341) experi¬ 
ence mean magnification > 1.4. Postage stamps of 
these dropouts are shown in Figure 6 and their lensing 
properties are presented in Table 2. As described in Sec¬ 
tion 2.1 the fiducial BoRG contamination fraction is 42% 
(Bradley et al. 2012; Schmidt et al. 2014a) meaning that 
some of the sources presented here may be lower redshift 
interlopers (e.g., Hayes et al. 2012; Bouwens et al. 2013). 
Without further photometry and/or spectroscopy we 
cannot identify which of the sources are interlopers, 
but we note that the photometric redshift PDFs for 
these four sources (obtained from BPZ) all have strong 
peaks at z ~ 8, suggesting a higher probability than the 
average (58%) for these particular objects to be true 
z ~ 8 sources. Interestingly, borg_1301±0000_160 is the 
brightest dropout in the survey, with mj = 25.5 ± 0.2, 
and appears tangentially elongated in the J-band image 
(middle panel of Figure 6). This object is also a very 
interesting target for further imaging and spectroscopic 
follow-up. 

We note that our method assigns a significantly 
lower velocity dispersion to the potential strong lens 
(borg_0440-5244_647) than the one estimated by Barone- 
Nugent et al. (2013) in their presentation of this object. 
They estimated the velocity dispersion of the deflector 
to be cr ^ 300 km s“^, whereas our method estimates a 
mean velocity dispersion of ^ 170 ± 33 km s“^. This is 
likely to be because our method does not account for lens¬ 
ing by groups and clusters, while Barone-Nugent et al. 
(2013) suggest that this dropout is lensed by a group of at 
least two objects at z ~ 1.8, of which borg_0440-5244_647 
is the largest. They estimated velocity dispersions of the 
deflector galaxies by using an abundance matching rela¬ 
tion between mass and luminosity, derived from Cooray 
(2005), and measuring the angular size of the lensing ob¬ 
jects. However, when using a redshift-dependent Faber 
& Jackson (1976) relation (Barone-Nugent et al. 2015) 
similar to ours (Table 1) they estimate the velocity dis¬ 
persion of this single galaxy to be ~ 180±46 km s“^ (via 
private communication), which agrees with our result. 
Neglecting group-scale lensing is a potential limitation of 
our method, which may underestimate magnification in a 
few cases. However the impact on the overall estimation 
of the LF inference is negligible since the phenomenon is 
so rare. 
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rriAB 

Fig. 5.— Evolving Faber-Jackson relation for massive galaxies with redshift. Data for z < 1 are from Treu et al. (2005) (red and green 
circles), data for z > 1 are from Belli et al. (2014a) (blue triangles) and Belli et al. (2014b) (blue crosses). Red points indicate apparent 
magnitude in the F606W band {z < 0.5), green points have magnitudes in the F098M band (0.5 < 2 < 1), and blue points are data with 
magnitudes in the F160W band {z > 2). Only galaxies with a > 200 km s“^ were used in the fitting. The slope of the relation in velocity 
dispersion and magnitude is fixed at the Faber Jackson (1976) result of L oc cr^. We fit the evolution with redshift, which changes the 
intercept of the line on the velocity dispersion axis. The uncertainty in magnitude is 0.1 mag which is a fiducial value given the fitting 
procedures. The black dashed lines shows the mean fit for the mean redshift of objects in each plotted bin. The fitting parameters are 
given in Table 1. 


TABLE 2 

Strong and intermediate tensing parameters derived by estimating velocity dispersions of bright foreground galaxies 

CLOSE TO 2 8 dropouts 


Field 

Dropout ID 

Jl25" 

Foreground ID 


Separation {") 

(Tinf(km s 

^BR (") 

IJ- 

borg_0436-5259 

1233^^ 

27.1 ± 0.2 

1191 

1.52 ± 0.03 

2.79 

294 ± 47 

1.32 ± 0.40 

2.05 ± 0.52 

borg_1301-|-0000 

160^ 

25.5 ±0.2 

144 

1.14± 0.15 

1.99 

184 ± 31 

0.60 ± 0.20 

1.47 ±0.30 

borg_1408-l-5503 

980'^ 

27.0 ±0.2 

959 

0.40 ± 0.06 

3.11 

193 ± 69 

1.01 ± 0.70 

1.54 ±0.62 

borg 2155-4411 

341'= 

26.6 ±0.2 

244 

0.74 ± 0.11 

2.27 

216 ±22 

0.97 ±0.20 

1.80 ±0.33 


Note. — ^ Total (AUTOMAG) apparent magnitude in the J-band of the dropout (Bradley et al. 2012). 
^ Strongly-lensed candidate. ^ 5cr source. ^ 8c7 source. 


5. WEAK LENSING tervening matter. We used the Pangloss code^ developed 

In this section we discuss the methods used to find the 1?^ Collett et al. ^013) that generates leasing parameters 
PDFs for magnification of a source at z ^ 8 by all in- reconstructed Imes-of-sight. We describe the produc- 


http://github.com/drphilmarshall/Pangloss 
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Fig. 6.— The four BoRG dropouts (from top left to bottom right: borg_0436-5259_1233, borg_1301+0000_160, borg_1408+5503 and 
borg_2155-4411_341) with significant magnification probabilities, shown in the F125W band with a Gaussian smoothing radius of 1 in 8” 
boxes. The solid red lines outline the dropouts with a 0.3” radius. The dashed green lines outline the potential foreground deflectors, with 
radius corresponding to the Einstein radius of an SIS deflector lensing a source at 2 = 8. The candidate strong lens system (borg_0436- 
5259-1233) is shown in the top left panel, and has an estimated magnification of /r = 2.05 it 0.52. Interestingly, borg_1301+0000_160 (top 
right) is the brightest dropout in the BoRG survey. The parameters for all of these objects are given in Table 2. 


tion of magnification PDFs from simulation data from 
the Millennium Simulation (Springel et al. 2005) in Sec¬ 
tion 5.1 and in Section 5.2 we present the BoRG field 
weak magnification PDFs. Our PDFs agree well with 
other theoretical work at lower redshifts (Hilbert et al. 
2007, 2009; Greene et al. 2013). 


5.1. Estimating Magnification from Simulation Catalogs 

The weak lensing reconstruction model developed by 
Gollett et al. (2013) takes simulation halo catalogs and 
places halos in a three-dimensional grid, with each halo 
contributing convergence Ki and shear 7 ^ along a line-of- 
sight to a source at a given redshift. Halos are modeled 
as truncated NFW profiles (Baltz et al. 2009); 


p(r) 


PNFw(r) 



(13) 


where we used the truncation radius rt = 5 r 2 oo, shown 
to be robust by Collett et al. (2013). Where r 2 oo is the 
radius at which the mass density falls to 200 times the 
critical mass density of the universe. The convergence 
and shear derived from this profile are given in Baltz 
et al. (2009). Magnification due to all intervening deflec¬ 
tors along a line-of-sight is given by Equation (7). 

We built PDFs for all lensing parameters by sampling 
over 10^ of lines-of-sight. As described in Section 2.3 
we used lightcones built from the Millennium Simulation 
(Henriques et al. 2012; Springel et al. 2005). 

The simulated catalogs provide a list of halos with as¬ 
sociated galaxies, but they do not include other dark 
structure, clumped in filaments and absent in voids. This 
missing matter will affect the overall density of the uni¬ 
verse so it is necessary to take this into account when 
estimating k and p,. We account for this by subtracting 
convergence from redshift slices so that the mean con¬ 
vergence along all lines-of-sight in the catalogs to a given 
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Fig. 7.— Comparison of the overdensity of lines-of-sight in the 
Millennium Simulation and the BoRG fields. ^ = Ui/ntot where 
is the number of objects per unit area above a certain flux limit in 
each lightcone (or real BoRG field) and ntot is the total number of 
objects above the same flux limit divided by the total survey area. 
We use a flux limit of m < 24 in F125W (J-band). 

redshift equals zero, and the mean magnification is unity, 
as they should be. 

Following work by Suyu et al. (2010) and Greene et al. 
(2013), we compare lines-of-sight in the BoRG fields with 
simulation data based on relative density of objects. We 
define the overdensity parameter 


where Ui is the number of objects per unit area in each 
lightcone (or real field) and ntot is the total number of 
objects divided by the total survey area. Given that 
the simulation catalogs are ~ 500 x larger than the total 
BoRG survey area we expect them to give representative 
results. 

We then calculate the number of objects per square 
arcsecond brighter than m = 24 in the J-band in each of 
the BoRG fields compared to the total number of objects 
above this flux limit in the whole survey. Similarly we 
calculate the overdensity of objects above the same limit 
in the simulated lightcones. Henriques et al. (2012) in¬ 
clude mock photometry based on stellar population syn¬ 
thesis codes by Maraston (2005) which include J-band 
magnitudes. As shown in Figure 7, the distribution of 
overdensities for the observed data is within the range 
of that for simulated data. Finally, to generate magni¬ 
fication PDFs for a given BoRG field, we combine the 
magnifications from all simulation lines-of-sight which 
are within ±2% in overdensity of the observed value. 

In Figure 8 we plot the magnification PDFs for a source 
at various redshifts over all lines-of-sight. As the source 
redshift increases, the peak of the distribution shifts to 
lower magnification, but the high-magnification tail be¬ 
comes more important, such that the mean magnification 
over all lines-of-sight remains unity. We match results for 
z < 6 from Hilbert et al. (2007) well. It is clear that there 
is little change in the distribution between zg = 6 and 
zs = 8, as there are negligible numbers of large halos 
above z > 5. 

In Figure 9 we plot the magnification PDFs for a va- 





Fig. 8.— Probability distribution function for magnification for 
four values of source redshift. The dashed line marks the mean 
magnification of the universe. These results compare well with 
Hilbert et al. (2007). Due to the lack of significant mass between 
2 ~ 6 and .2 ~ 8 there is little change in the distributions of mag¬ 
nification for sources at those redshifts, as the total convergence 
does not change much. 





Fig. 9.— Probability distribution function for magnification for 
a range of values of overdensities for a source at 2 = 8. More over- 
dense lines-of-sight are skewed towards higher magnification, with 
a broad distribution. More underdense lines-of-sight are skewed 
towards lower magnification, with a narrower distribution due to 
the deficit of intervening mass. 

riety of overdensities. The more overdense lines-of-sight 
produce a higher mean magnification, as expected, but 
also have a greater variance than the distributions for 
underdense lines-of-sight. This agrees well with the esti¬ 
mates at lower redshift by Greene et al. (2013). 

5.2. BoRG Weak Lensing Magnification PDFs 

The kernel density estimates (Rosenblatt 1956; Parzen 
1962) fit to the magnification PDFs for all the BoRG 
fields are shown in Figure 10. As expected, the BoRG 
fields do not have significant over- or underdensities, but 
are rather typical of blank fields at z ^ 8, as shown in 
Figure 8. 

There is significant motivation for the magnification 
PDFs to take a log-normal form. The 3D matter den- 
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Fig. 10. — Probability distribution function for magnification for 
all of the BoRG fields, with a source at 2 ~ 8. The lines are kernel 
density estimations to the distributions. It is clear there is little 
range in overdensity for the BoRG fields. 


sity distribution of the universe is well-described by a 
log-normal random field (Coles & Jones 1991), and weak 
leasing probability distributions arise directly from the 
mass distribution. However, when accounting for the 
magnification bias in individual fields to infer the LF 
from the dropout sample (see Section 6) it was necessary 
to express the magnification distributions in a form that 
could easily convolve analytically with a Gaussian dis¬ 
tribution (for more details see Appendix A). For this we 
used a Bayesian MCMC approach to fit the distributions 
of magnification for each field as a linear sum of Gaussian 
functions. 


6. RECOMPUTING THE LF 

In this section we outline the method of estimating 
the z ~ 8 LF from the BoRG high-redshift candidates, 
taking the magnification bias into account. 

Following Schmidt et al. (2014a), who did not account 
for the magnification bias when estimating the BoRG 
z ~ 8 LF, we use the Bayesian inference method devised 
by Kelly et al. (2008), which is described in Section 6.1 
and in Appendix A. In Section 6.2 we describe in more 
detail how we take into account the weak and intermedi¬ 
ate leasing magnification. 

6.1. Bayesian Estimation of the LF 

As in Schmidt et al. (2014a), we assume that the in¬ 
trinsic luminosity function is modeled by the Schechter 
function in Equation (3). In order to facilitate compari¬ 
son with our previous work we use the sample of 38 BoRG 
Y-band dropouts and 59 additional fainter dropouts from 
the Hubble Ultra-Deep Field (HUDF) and Early Release 
Science (ERS) programs (Bouwens et al. 2011). 

Bayesian statistics allows us to express the posterior 
probability that the LF is fit by a Schechter function 
with parameters 9 = (a, L*,4'*) given the observed lu¬ 
minosity Lj,obs of the dropouts in the J-band, and the 
non-detections in the V-band (ly = 0), as the product of 
the prior on the Schechter parameters and the likelihood: 

p{0 I Lj.obs, ly = 0) oc p{9) X p(Lj,obs, /y = 0 I 0) (15) 


This posterior probability can be expressed (see Ap¬ 
pendix A for full details) as: 


p{o I Lj,obs, /v = 0) cx p{6) X X c;, 

C r 


rN^ 


n 


1 - 


A, 


Pl-^sky 
n 

YIp^Lj ,obs,2 \d) 


p{I=l\ 9) 




(16) 


Where we iterate over I fields with i z ~ 8 candidates. 
Here is the number of high-z dropouts in the sur¬ 
veyed comoving cosmological volume, Ai is the area of 
the individual C fields in Schmidt et al. (2014a), which 

each contain ci high redshift candidates (n = c;). 

Each candidate has an assumed contamination fraction 
of fi- We use a fiducial value for the contamination of 
42% for the BoRG sample, the contamination fractions 
for the HUDF/ERS samples are included in the selec¬ 
tion function, see Appendix A, as described in Oesch 
et al. (2012); Bradley et al. (2012) and Schmidt et al. 
(2014a). Ghanging the contamination value in the range 
/ = 0 — 0.60 effects the characteristic magnitude and the 
number density of the LF by less than their estimated 
1(7 uncertainties, and the change in the faint-end slope 
is comparable to its la uncertainty (Bradley et al. 2012; 
Schmidt et al. 2014a). The Bayesian framework allows 
us to accurately estimate the LF parameters accounting 
for contamination. Agky is the area of the full sky The 
factors are binomial coefficients which are the fully 
correct method of modeling source counts. 

We assume uniform priors on a, log^Q L* and log^Q N^. 
p{I = 11 6>) is the probability distribution of an object 
making it into the dropout sample based on the pho¬ 
tometric selection described in Schmidt et al. (2014a). 
p(Lj_obs.i I 9) is the likelihood function for the observed 
J-band luminosity of the Fth object in the sample. 

The last term includes marginalization over the mag¬ 
nification PDF: 


P(Aj,obs I 9') — J' j p{p) p(Lj^obs I /^Lj true) 

^ P(Aj^true I ^) true (1^) 


In Appendix A we give the expanded expression of the 
posterior distribution from Equation (16) used when per¬ 
forming the LF parameter inference and describe the 
derivation and motivation for Equation (17). We refer 
to Appendix A and Schmidt et al. (2014a) for further 
details. 


6.2. Including the Lensing Corrections 
6.2.1. Analytic Form for Magnification PDFs 

In order to make integration of Equation (17) compu¬ 
tationally feasible we require a simple analytic form for 
p{p) that will convolve simply with a Gaussian distri¬ 
bution (see Appendix A). As described in Section 5.2, 
the weak lensing magnification PDF is well-fit by a log¬ 
normal distribution. However, this cannot be convolved 
analytically with a Gaussian. 

Therefore, we fit the magnification PDFs from all 
regimes as a linear combination of Gaussian functions 
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with different means and standard deviations. The weak 
leasing magnification PDFs (see Section 5.2) are well-fit 
by a combination of three Gaussian functions. The inter¬ 
mediate leasing PDFs (see Section 4.2) are also well-fit 
by a combination of three Gaussian functions. 

6.2.2. Combining Lensing Regimes 

All of the fields have a weak lensing magnification 
PDFs based on their overdensity (see Section 5), but we 
have also identified one strongly-lensed candidate and 
three dropouts close to large foreground galaxies that 
produce an intermediate magnification PDF (see Sec¬ 
tion 4.2). 

To account for the magnification bias, we need to use 
the correct magnification PDF for each field. In the case 
when a strong or intermediate lens appears present, we 
split the field into two parts for the calculation of the 
posterior: one is a circle with radius 10 6*er containing 
the dropout and the deflector, where we use the strong or 
intermediate lens magnification PDF. For the remainder 
of the field we use the weak lensing magnification PDF. 

Whilst the total flux across the sky is conserved, lo¬ 
cally over- or underdensities that produce magnification 
not only magnify fluxes, but also increases areas. Hence, 
the individual BoRG fields we observe have been magni¬ 
fied (or demagnified) from their true sizes. We account 
for this in the posterior probability Equation (16) by di¬ 
viding the measured area of each field by the mean mag¬ 
nification in that field, /I; from the magnification PDFs. 
For weak lensing magnification PDFs (/r;) ^ 1. For the 
intermediate lensing case 1.4 < <2 due to our selec¬ 

tion process. 

As magnification is most important for the bright-end 
of the LF, and negligible at the faint end, for simplicity 
and without loss of precision, we adopt /x = 1 for the 59 
fainter dropouts (Bouwens et al. 2011). Additionally, one 
of the BoRG fields (borg_1815_3244) is centered on the 
Galactic plane and is dominated by stars. We discard 
this field in our calculation of the LF. 

7. RESULTS 

Using the framework described in Section 6 to account 
for the magnification bias we present our estimation of 
the z ~ 8 galaxy LF based on our sample of 97 z ~ 8 
LBGs (described in Section 2). First, in Section 7.1 we 
compare our estimates of strong and intermediate lens¬ 
ing probabilities with the actual observations. Then, in 
Section 7.2 we carry out the inference of the z ~ 8 LF. 
Finally, in Section 7.3, we use our semi-analytical model 
of strong lensing optical depths described in Section 4 to 
predict the form of observed LFs at z > 8. 

7.1. Strong and Intermediate Lensing Events in the 
BoRG Survey 

The simple SIS strong lensing model described in Sec¬ 
tion 4.1 predicts the probability of z ~ 8 sources in the 
BoRG survey being multiply imaged to be ~ 3—15%, in¬ 
creasing as the field limiting magnitude becomes brighter 
than M*. The majority of the BoRG fields have a 
multiple-image probability for high-redshift sources of 
< 10% (see Figure 4). We predict that 1-2 of the 38 
BoRG Y-band dropouts may be strongly lensed. 

One candidate strong lens system in BoRG was pre¬ 
sented by Barone-Nugent et al. (2013), a rigorous search 


for strong lenses in all 71 BoRG fields as part of this work 
revealed one more candidate. Additionally, this search 
revealed three candidate intermediate lens systems, with 
fj, > 1.4. These candidates are presented in Figure 6 and 
Table 2. Whilst strong lensing creates larger magnifica¬ 
tion, the probability of encountering a strong lens along 
the line-of-sight is low: as shown in Figure 3 the optical 
depth is roughly r « 0.31% for a source at z = 8. The 
optical depth for intermediate lensing is much higher: for 
an object to experience intermediate lensing it must be 
within 3.50ER of the foreground deflector, resulting in 
r ss 4% for a source at z = 8. Thus, intermediate lens¬ 
ing offers an additional boost to the flux of high-redshift 
galaxies, and must be correctly accounted for in estima¬ 
tions of the LF. 

7.2. Inference of the Intrinsic z ~ 8 LF 

We estimate the z ~ 8 LF from the sample of 97 LBG 
described in Section 2, including the 38 S/Nj > 5 objects 
from the BoRG survey, including the effects of magnifica¬ 
tion bias. We sample the posterior distribution function 
for the Schechter function parameters with an MGMG 
chain of 40 000 steps. 

The results of the estimated LF are shown in Figure 11, 
and the correlations between the Schechter function pa¬ 
rameters and their PDFs are shown in Figure 12. We 
plot the results of Schmidt et al. (2014a) for comparison 
in both figures. We see a small deviation from the un¬ 
corrected LF of ~ 0.15 mag at the limit of the brightest 
BoRG source, and there is negligible difference between 
the LFs at M > —21. The Schechter function parame¬ 
ters for the new LF are within the uncertainties of the 
estimation by Schmidt et al. (2014a), though we find a 
slightly fainter value of M* and higher value of 4** than 
Schmidt et al. (2014a). This is expected because of the 
slight deviation at the bright end of the LF, and there is 
a strong correlation between these parameters, as shown 
in Figure 12. It is clear that magnification bias is not 
a significant effect at this redshift and the luminosity 
range of the BoRG sources. This also demonstrates that 
although we predict 3 — 15% of the BoRG sources are 
strongly lensed this does not affect the LF within the 
survey limits, as predicted in Section 3.4. 

Our results are in good agreement with those of Fialkov 
& Loeb (2015) who use an independent semi-analytic 
method to show the effect of magnification bias is small 
below M > —21.5. Fialkov & Loeb (2015) predict that 
if the brightest observed galaxy has absolute magnitude 
Muv = —24.5 - lying in the significantly distorted tail 
of the magnified LF (Wyithe et al. 2011) - there is a 
~ 13.3% discrepancy in the normalization of a Schechter 
LF at z ~ 8 for sampled galaxies with /Xmax = 2 (i.e. only 
weak and intermediate lensing effects) compared to the 
intrinsic LF. Whilst this upper limit is several orders of 
magnitude brighter than currently observed, this demon¬ 
strates that it will be important to include the effects of 
magnification bias from weak and intermediate lensing 
in surveys that find extremely bright galaxies. 

Table 3 summarizes the estimated Schechter function 
parameters for this LF in comparison with other recent 
LF estimates from the literature. We find that our fit 
parameters are in good agreement with the recent liter¬ 
ature, demonstrating that magnification bias is not af¬ 
fecting current z ~ 8 LF observations. Note that our 
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Muv 

Fig. 11.— The intrinsic 2: ~ 8 LF, which is well-described by 
a Schechter (1976) function, including the magnification bias due 
to weak and intermediate lensing in all BoRG fields (solid black 
line). We plot the LF without the treatment of the magnifica¬ 
tion bias (Schmidt et al. 2014a) for comparison (dashed red line). 
The lines corresponds to the median values of the MCMC samples 
and the shaded regions correspond to the 68% confidence region of 
the samples. The LF estimated here is virtually indistinguishable 
from that of Schmidt et al. (2014a), demonstrating that magni¬ 
fication bias is not a significant effect at 2 ~ 8. The Schechter 
parameters for this LF are given in Table 3 along with literature 
values. The binned data from BoRG12 (Bradley et al. 2012) and 
the faint HUDF/ERS candidates (Bouwens et al. 2011) are also 
plotted as blue and green points respectively. The inverted green 
triangle denotes the brightest BoRG dropout. We note that the 
LF is estimated from the unbinned data. 

results have significantly smaller error bars than those of 
Finkelstein et al. (2014), because their sample contains 
only 3 z ^ 8 galaxies brighter than M = —21, making 
their fit less well-constrained at the bright end. 

Our results show that magnification bias does not af¬ 
fect current estimates of the LF at z < 8 and therefore 
cannot explain the apparent flattening of the bright-end 
of the LF recently observed by Bowler et al. (2014a,b) 
and Finkelstein et al. (2014) at z ^ 7 — 8. Bowler 
et al. (2014a,b) accounted for strong lensing of bright 
sources, but they still find a deviation of ~ 0.4 mag from 
a Schechter fit at M = —22. We predict a lensed frac¬ 
tion of ^ 3 — 15% for bright galaxies (Figure 4) from the 
BoRG survey which is essentially free of cosmic variance 
(Trenti & Stiavelli 2008), so providing cosmic variance 
and contamination by lower redshift interlopers (Hayes 
et al. 2012; Bradley et al. 2012; Bouwens et al. 2013; 
Schmidt et al. 2014a) were correctly accounted for in the 
work of Bowler et al. (2014a,b) and Finkelstein et al. 
(2014), we expect the magnification bias to be negligible 
in the bright-end of these LFs. This lends credence to the 
interpretation that these observations may be the result 
of the changing intrinsic properties of galaxies at z > 7, 
possibly due to changing dust fractions (Cai et al. 2014) 
and/or feedback processes (Somerville et al. 2012). 

7.3. Predictions for z > 8 and Future Surveys 

There is clear evolution in the LF for z < 8 (e.g. 
Bouwens et al. 2007; van der Burg et al. 2010; Bouwens 
et al. 2014; Bowler et al. 2014a; Finkelstein et al. 2014), 
and this is expected to continue to higher redshifts. How¬ 
ever, the processes which drive this evolution are not 


well-understood: the evolution is thought to follow hi¬ 
erarchical structure formation and the evolution of the 
halo mass function (Vale & Ostriker 2004), but there 
are also important quenching processes that may reduce 
star formation in massive galaxies (Schneider et al. 2006; 
Somerville et al. 2012), and changes in the amount of 
dust present in galaxies will affect the attenuation of flux. 
Thus there are a multitude of theoretical models for the 
evolution of the LF. 

The gravitationally lensed LF (Equation (10)) exhibits 
a significant ‘kick’ in the bright-end tail for M < — 22 at 
z ^ 8. This is just beyond the brightest BoRG objects, so 
it is unlikely that the BoRG survey observes the regime 
of magnification bias at the bright-end. This is in agree¬ 
ment with theoretical studies by Wyithe et al. (2011) and 
Fialkov & Loeb (2015). However, in upcoming wide-area 
surveys magnification bias presents a useful tool to test 
LF evolution models because it allows us to probe the 
bright end, where there are large theoretical uncertain¬ 
ties and the evolution is expected to be fast (Bowler et al. 
2014a). 

In order to explore the range of possible scenarios, in 
Figure 13 we plot the predicted intrinsic (dashed lines) 
and observed (solid lines) LFs for a range of redshifts, 
comparing a variety of evolution models. We assume 
these models are the intrinsic LFs at a given redshift 
and used Equation (10) to estimate the observed LE. We 
plot the BoRG z ~ 8 LF (Schmidt et al. 2014a) for com¬ 
parison. Additionally, we mark the comoving volumes 
and magnitude ranges accessible to future high-redshift 
surveys. 

The top left panel shows the LF model from Bouwens 
et al. (2014) which is an extrapolation from observations 
at z < 10. The top right panel shows the LF model from 
Finkelstein et al. (2014) which is an extrapolation from 
observations at 4 < z < 8. The bottom left panel shows 
the model developed by Munoz (2012) which follows the 
evolution of the halo mass function, and includes dust 
attenuation. The bottom right panel is a model from 
Behroozi & Silk (2015) constructed from a comparison of 
the specific star formation rate to the specific halo mass 
accretion rate, and including dust models from Gharlot 
& Fall (2000). The four models have significantly dif¬ 
ferent behaviors at the bright end. While the Bouwens 
et al. (2014) model has by construction a bright end that 
is very similar to that measured at lower redshifts, the 
Munoz (2012) model has a very shallow bright end, and 
the Finkelstein et al. (2014) and Behroozi & Silk (2015) 
models are in-between. As a result, the effects of mag¬ 
nification bias (which are stronger for the steeper LF) 
are very different: negligible in the Munoz (2012) case 
and appreciable in the three other cases. However, the 
bright end of the Munoz (2012) model is the easier one 
to test observationally, within reach of a James Webb 
Space Telescope medium depth, medium width survey 
(e.g. JWST MD Windhorst et al. 2006). 

Except in the case of a very shallow bright end, we 
do not expect the magnification bias to be significant in 
our upcoming BoRG z ~ 9,10 survey (HST Gycle 22, 
PI Trenti). In all cases, it is clear that surveys covering 
> 100 deg^, e.g. Euclid and WFIRST, should find many 
bright z > 8 LBGs. We expect the observed high-redshift 
galaxy samples will be dominated by magnification bias 
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Fig. 12.— The correlations between the 2 ~ 8 LF Schechter function parameters (a, M* and ^*) estimated from the BoRG dropouts 
including treatment of magnification bias (black), compared to the parameters obtained without the treatment of magnification bias (red, 
Schmidt et al. 2014a) with \g and 2(7 confidence contours. There is clear correlation between all three parameters. The top panels show 
the marginalized PDFs for each parameter. 


TABLE 3 

Comparison of 2 : ~ 8 Schechter LF Parameters 


Reference 

M* 

a 

logio'J'* [Mpc 

This work 

Finkelstein et al. (2014) 
Bouwens et al. (2014) 
Schmidt et al. (2014a) ba 
Schmidt et al. (2014a) 8 (t 
M cLure et al. (2013) 
Schenker et al. (2013) 

-19.85:^;^^ 
-20.89^;^^ 
-20.63 ±0.36 
-20.151°^® 
—20 40"^®'^^ 

zu.^u_o.55 

—20 

ZU.iZ_o 48 
zu.^4 0.35 

-^•'^-0.29 

—2 

^•■^“-0.40 
-2.02 ± 0.23 

1 o'7+0.26 

-0.26 
^•0o_Q 29 

9 n9+0.22 

'^•'^^-0.23 

-^•^^-0.24 

— 3 00"^*^'^^ 

a.uu_o 31 

—4 14+0-®5 

-3.68 ±0.32 

^•^^-0.34 

3 K-| +0.36 

3 

o.oo_Q 47 

a.ou o 32 


in these surveys. We predict almost all z ~ 8 sources 
in Euclid will have been strongly lensed. The framework 
developed in this work will be crucial for determining the 
intrinsic luminosity of high-redshift sources found in such 
surveys. 

Our results confirm the suggestion by Wyithe et al. 
(2011) that magnification bias will be important to probe 
the bright end of the LF at high redshift. However, we 
find that the magnitude of the effect is less pronounced 
than in that study, owing mostly to our accounting for 
the redshift evolution of the deflector population. 


8. SUMMARY AND CONCLUSION 

We have introduced a systematic way to account for 
the magnification bias in estimations of high-redshift 
LFs. The method involves estimating the probability 
density function for weak leasing magnification along a 
given line-of-sight by comparison with results from the 
reconstruction of simulated halo data, and by estimating 
the strong and intermediate leasing magnification PDF 
of dropouts due to massive deflector galaxies in close 
proximity to the dropout. 
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Fig. 13.— Predicted observed LFs for z > 8 redshifts. For 2 = 8 we use the Schechter LF from Schmidt et al. (2014a), plotted as a 
thick black line. The white band indicates the error on the Schechter function parameters, and the thin black line is the extrapolation of 
the LF beyond the observational limit. We show the regions of magnitude and volume observable by current and future surveys: the total 
BoRG survey including the 2 ~ 8 survey described in Section 2.1 and the upcoming BoRG 2 ~ 9,10 survey (HST Cycle 22, PI: Trenti); the 
James Webb Telescope Medium Deep (JWST MD) (Windhorst et al. 2006); the Wide-Field Imaging Surveyor for High-Redshift Ultra-Deep 
Field (WISH UDF, http://wishmission.org/en/doc.html); the Wide-Field Infrared Survey Telescope High Latitude Survey (WFIRST HLS) 
(Spergel et al. 2013) and the Euclid Wide Survey (WS) (Laureijs et al. 2011). As explained in the text, BoRG does not survey enough area 
to observe the rarest bright sources which are most affected by magnification bias, but future wide-held surveys will be dominated by this 
effect. (Top Left) For 2 > 8 we use the LF model from Bouwens et al. (2014) which is an extrapolation from 2 ~ 10. (Top Right) For 
2 > 8 we use extrapolate the evolution of the Schechter function parameters over 4 < 2 < 8 from Finkelstein et al. (2014). (Bottom Left) 
For 2 > 8 we use the luminosity model from Munoz (2012) which is based on the evolution of the halo mass function. These do not exhibit 
the sharp cut-off at the bright-end and are not affected by magnihcation bias. (Bottom Right) For 2 > 8 we use the LF evolution model 
from Behroozi Sz Silk (2015). The dashed lines indicate the intrinsic LFs, the solid lines are the observed LFs including the magnification 
bias calculated using Equation (10). 
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We applied this method to estimate the 2 ~ 8 LF from 
the 38 BoRG Y-band dropouts and 59 fainter dropouts 
from Bouwens et al. (2011). Our main results are sum¬ 
marized as follows: 

(a) The probability of a BoRG z ~ 8 dropout being mul¬ 
tiply imaged is ~ 3 — 15%, increasing with limit¬ 
ing magnitude. This is consistent with finding two 
strongly-lensed dropouts in the BORG survey: the 
candidate system presented in Barone-Nugent et al. 

(2013), and the additional strongly-lensed candidate 
dropout in this paper. We also find three dropouts 
which may experience significant magnification with¬ 
out multiple imaging, consistent with our expecta¬ 
tions. 

(b) We extended the Bayesian formalism for the estima¬ 
tion of the LF parameters presented by Schmidt et al. 

(2014a) to account for the magnification bias. This 
involves marginalizing over the magnification PDFs 
for strong and weak lensing effects. The inferred 
Schechter function parameters are: 

= -19.85t°jo, 

O' — —1 72+^-^*^ 
a — i. < ^_o. 29 ’ 

logio 'I'*[Mpc-3] =-3.00™, 

These values do not differ significantly from estimates 
not accounting for the magnification bias. 

(c) Thus magnification bias cannot be an explanation 
for the apparent flattening of the bright-end of the 
LF recently observed by Bowler et al. (2014a,b) and 
Finkelstein et al. (2014). 

APPENDIX 

A. BAYESIAN FRAMEWORK FOR ESTIMATING THE LUMINOSITY FUNCTION 

We use Bayesian statistics to find the relationship between the prior probability of the z ~ 8 dropouts being galaxies 
with LF Schechter parameters 6 = (a, L*, dt*), and these parameters’ posterior probability given the dropout 
candidates’ detection threshold in the J-band, assuming their non-detection in the V-band. The posterior probability 
is given by: 

p{0 I Lj,obs, /y = 0) oc p{0) X p(Lj,obs, Ly = 0 I d) (Al) 

where the last term is the likelihood and p{9) is the prior on the LF parameters. We will assume uniform priors on a 
and logiQ L*. 

We can expand the expression for the posterior: 

C riz 

p{0 I Tj.obs, /v = 0) (X p{9) J] [1 - Ai/A,^y p{I = 110)]"^^-='^ x J]p(Lj,obs,* | 9) 



X n [1 - p{I = l\9)f^-‘^‘^ X nP(^J.obs.* I 9) (A2) 

I i 

where the terms are binomial coefficients which correctly model the distribution of source counts. and Nc are 
the number of high-redshift sources given the intrinsic LF and the number of potential contaminants in the Universe 
respectively. We will assume a uniform prior on logj^Q Nz- In the observed sample the number of high-redshift sources 
and contaminants are given by and Uc- The total number of galaxies in the observed sample, rit is given by their 
sum. We take the product over C individual observed fields where c/ represents the number of galaxies in the Fth field 
with rit also given by the sum of c/ over all of the fields. The fraction of the sky covered by the Fth field is given by 
Ai/Asky. The contamination fraction in each field, fi is set at the fiducial value of 42% (Schmidt et al. 2014a; Bradley 
et al. 2012) for the BoRG sources, the contamination fraction for the fainter HUDF/ERS sources (Bouwens et al. 2011) 
is included in the selection function (see below). 


(d) The z ~ 8 LF appears significantly magnified for ex¬ 
tremely bright galaxies (Muv < —22). Though cur¬ 
rent surveys have not observed such rare, luminous 
galaxies, future wide-held surveys will probe this re¬ 
gion. For surveys > 100 deg^, e.g. WFIRST, Euclid, 
we predict that samples of z > 8 galaxies will be 
dominated by magnihcation bias. 

(e) Magnihcation bias will be a useful tool to distinguish 
between high-redshift LF evolution models. In par¬ 
ticular it could help determine whether the LF tran¬ 
sitions from a Schechter form to a power-law form 
at high redshift, indicating signihcant changes in the 
astrophysical properties of those galaxies. 
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The last term in Equation (A2) is the likelihood for the i’th object in the sample. In Schmidt et al. (2014a) this was 
expressed as: 

pOO 

P(-^J,obs I ~ / P(-^J,obs I -^J,true) p{^J,true \ ^-^J,true (-^^) 

^0 

•^(-^J,obs I -^J,true5 field) | -L ) <i-£/J^true 

where we use p{L \ 6) cx (see Equation (1) of Kelly et al. (2008)). The function gamma(Lj^true | a^L*) is related 

to the Schechter LF (Equation (3)) as gamma(L | a, L*) = • 



(Tj^obs I Tj^true 1 geld) — 


dTj geldC^TT 


exp 


(l^J,obs l^J,true) 


2 <5^1: 


field 


(A4) 


represents the true luminosity inferred from the observations assuming a Gaussian measurement error with dAj,field 
being the median photometric error in the J-band in the given field. 

In order to include the effects of the magnification bias, we must integrate over the nuisance parameter Aj,mag) which 
represents the luminosity of an object in the J-band, magnified above its true luminosity. Including this Equation (A3) 
becomes: 


p(Aj.obs \0)=J J p(Lj,obs I O ,mag ) p(o ,mag lO jtrue ) PiLj drue I 9) dLj drue 


(A5) 


where p(Aj,obs | J^j,mag) is now the term with Gaussian measurement errors similar to Equation (A4), given that we 
make observations of magnified luminosities: 


(Aj^obs I Tj^magi geld) — 




exp 


(^J ,obs Tj pnag) 

2 ggij 


(A6) 


To find the probability that luminosity is magnified from its true luminosity, p(Aj,mag I Aj^true)) we must integrate 
over the full magnification probability density: 


p(Aj ,mag |Aj jtrue) — ,mag \P, Aj ,true ) P{p) dfi 

We can marginalize over Aj^mag in the first part of Equation (A5): 

l^(Aj,obs I Aj true) — J' 4*(Aj^obs I Aj uiag) P{d^J^niag \ Aj^true) dAj ^lag 

= J j P(^^ ,mag Im, Aj ,true) pi.p^ p(Aj^obs 1 Aj juag; l^Aj geld) df-t dAj rng 
— J' J' '^(Aj^niag /^Aj true) P{p^ A^(Aj obs | Aj rnag, l^Aj geld) d/i dAj , 
— J' P{p) A^(Aj obs I MAj true, dAj geld) dfi 
_ J‘ ^ f (Aj,obs /^Ajtrue) 


(A7) 


dAj^eldV^ 


2 geld 


dfj, 


(AS) 


Here we have used the Dirac delta function d(Aj_mag~MAj,true) to map true luminosities to magnified luminosities. To 
make computation of Equation (A2) feasible, we integrate Equation (A8) analytically and want to remove any Aj mag 
dependence we fit the magnification PDFs as a normalized linear combination of Gaussian terms with coefficients j3i 
centered on Jt, mag with standard deviation cri^mag : 


pip-) = 


1 


^ O' i,magA/ ^TT 

Equation (A8) can then be integrated analytically: 

|d^ 1 


exp 


(M - A*i,mag)^ 
2'^i%ag 


P(Aj,obs I Aj true) — ^ ^ 


ri I xt2 

Y ^i,mag-^J,true ' 


: exp 


field 


(Tj,obs Mi,mag'^J drue) 

2 (f^fimagAyrue + 


J,field 


(A9) 


(AIO) 
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As solid angle is also magnified in gravitational lensing we must divide the measured field area A; by the average 
magnification in each field JLi- If > 1 the fields we observe appear larger than their true sizes. 

We can therefore express Equation (A2) as (see Schmidt et al. (2014a) for details): 


p{0 I Aj.obs, Iv = 0) (X p (0) X 


n 


1 - 


A, 


r<X) /*oo 




/O ^0 


^(-^J,obs,l) ‘^(-^J,obs,/5-^J,true,/) 


1 


X IT / -^(Aj.obs.i, Aj ,trueT ) dLj ,true,i 

. Jo 


(All) 


Here we have defined A'(Lj_obs, Aj true) = p(Aj^obs I Aj^true) gamma(Lj true I ce, L*) and included the selection function 
5(Aj obs)- The selection function estimates the completeness of the source selection and has been obtain for each 
individual BoRG field as explained in Oesch et al. (2012); Bradley et al. (2012) and Schmidt et al. (2014a). 

Thus, Equation (All) is the posterior probability distribution for a sample of rit binomially distributed objects, 
assumed have an intrinsic Schechter LE of the form shown in Equation (3). The observed luminosity of each object is 
related to its true luminosity via a magnification PDF and an assumed Gaussian error distribution. 
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